Nonlinear screening of charge impurities in graphene 
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It is shown that a "vacuum polarization" induced by Coulomb potential in graphene leads to 
a strong suppression of electric charges even for undoped case (no charge carriers). A standard 
linear response theory is therefore not applicable to describe the screening of charge impurities in 
graphene. In particular, it overestimates essentially the contributions of charge impurities into the 
resistivity of graphene. 
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Graphene is a name given to an atomic layer of car- 
bon atoms packed into a hexagonal two-dimensional lat- 
tice. This term is widely used to describe the crystal 
structure and properties of graphite (which consists from 
graphene layers relatively loosely stacked on top of each 
other), carbon nanutubes, and large fullerenes. Very re- 
cently, a way has been found to prepare free-standing 
graphene [H, that is, real two-dimensional crystal (in 
contrast with numerous gwasj-two-dimensional systems 
known before). The graphene turns out to be a gap- 
less semiconductor with a very high electron mobility 
which makes it a perspective material, e.g., for ballistic 
field-effect transistor 0- It has been shownjj, 0] that the 
charge carriers in graphene are massless Dirac fermions 
with effective "velocity of light" of order of 10 6 ms . 
Due to this unusual electronic structure graphene demon- 
strates exotic transport properties, such as a new kind of 
the integer quantum Hall effect with half-integer quanti- 
zation of the Hall conductivity [H, 0, 0, 0, 0, 01, or finite 
conductivity in the limit of zero charge-carrier concen- 
tration H la EE i- 

One of peculiar transport properties of graphene, a mo- 
bility which is almost independent on the c harg e carrier 



tie cnargc 

3d Ha « 



concentration[3(, was explained in Refslia, lia as a re- 
sult of electron scattering by charge impurities. How- 
ever, a linear-response theory was used to take into ac- 
count screening effects. Rigorously speaking, this theory 
can be applied only assuming that the impurity potential 
is small in comparison with the Fermi energy; however, 
even in semiconductors where this condition can be, in 
general, broken this theory can be normally used and 
gives reasonable results (for the case of two-dimensional 
electron gas, see for review Refll4T). 

In this paper we calculate nonlinear screening of charge 
impurities in graphene taking into account a "vacuum 
polarization" effect in a region of strong potential. 

A general nonlinear theory of screening in the system 
of interacting particles can be formulated in a framework 
of the density functional approach 15J. In this theory a 
total potential V (r) acting on electrons equals 

V (r) = V (r) + V ind (r) (1) 

where Vq (r) is an external potential and Vind (r) is a 



potential induced by a redistribution of electron density: 
„2 r „ 



Vmd (r) = - / dv' n ^' ] " + V xc (r) . 
e J r-r' 



(2) 



where the first term is the Hartree potential and the sec- 
ond one is the exchange-correlation potential. We will 
consider here explicitly only a redistribution of charge 
carriers in the external impurity potential 
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(3) 



taking into account contributions of crystal lattice po- 
tential and of electrons in completely filled bands via di- 
electric constant e and compensated homogeneous charge 
density —en; for the case of graphene on quartz one 
should choose (l3j| e ~ 2.4 — 2.5. Here Z is the dimension- 
less impurity charge (to be specific, we will assume Z > 0; 
it can be easily demonstrated that, actually, in our final 
expressions Z should be just replaced by \Z\). This kind 
of approach is valid at a space scale much larger than a 
lattice constant; in all other aspects, it is formally exact 
until we specify the expressions for V xc and n [V (r)]. 
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A dimensionless coupling constant a = (where 
vf — 10 6 ms _1 is the Fermi velocity in graphene) deter- 
mining the strength of interelectron interactions is of or- 
der of 1 which means that it is probably hopeless to con- 
sider the many-particle problem for graphene quite rig- 
orously. We will use the Thomas- Fermi theory [10] which 
is, actually, the simplest approximation in the density 
functional approach. It is based on the two assumptions: 
(i) we neglect the exchange-correlation potential in com- 
parison with the Hartree potential in Eq.((2|) and (ii) we 
put n (r') = n [fi — V (r')], n (//) being a density of ho- 
mogeneous electron gas with chemical potential //. The 
former assumption means that we are interested in the 
long- wavelength response of the electron system and thus 
the long-range Coulomb forces dominate over the local 
exchange-correlation effects. The later one holds pro- 
vided that the external potential is smooth enough. A 
rigorous statement is that an addition of constant poten- 
tial is equivalent to the shift of the chemical potential. 
In particular, the Thomas-Fermi theory gives an exact 
expression for static inhomogeneous dielectric function 
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e(q) in the limit of small wavevectors q — > 0(171 The 
Thomas-Fermi theory of atoms is asymptotically exact 
in the limit of infinite nuclear charge 16] . Here we will 
use it just for semiquantitative analysis of the problem. 
In the Thomas- Fermi theory Eq.@ reads 



V ind (r) = - 



ir/ n[ M -F(r')]-nM (4) 



The function n (fi) is expressed via the density of states 
N(E): 

n(fi) = J dEf (E) N{E) = J dEN (E) (5) 

where / (E) is the Fermi function and the last equality is 
valid for zero temperature (we will restrict ourselves here 
only by this case). For the case of graphene with linear 
energy spectrum near the crossing points K and K' one 
has 
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where we have taken into account a factor 4 due to two 
valleys and two spin projections. 

Let us start first with the case of zero doping (fi = 0) 
where, according to the linear response theory, there is 
no screening at all. Substituting Eqs.fi]), ([3]), and ([6|) 
into Eq.(fT]), introducing the notation 



V(r) = —F(r) 
er 



(7) 



and integrating over the polar angle of vector r', we ob- 
tain a nonlinear integral equation for the function F (r) : 



F(r) - Z -*f*L' K ^ W> (8) 
7r J r r + r \ r + r I 
o v ' 



where 
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K{k) 



dip 
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(9) 



sin if 



is the complete elliptic integral, 



Q = 2 
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(10) 



for the case of graphene on SiC>2 Q — 2. 

We will see below that, actually, the integral in right- 
hand side of Eq.® is divergent at r' — 0; the reason 
is that the expression |6]) with the replacement of /x by 
V (r) is not applicable for a very small distances when 
the potential becomes comparable with the conduction 
bandwidth; thus we should introduce a cut-off at r' ~ a 
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FIG. 1: Graph of the function <j>(t) (Eg. fTfy ) 



where a is of order of a lattice constant. An exact value 
of a is not relevant, with a logarithmic accuracy. 

To proceed further we replace variables in Eq.([5]), 
r' = re', and introduce a notation F (lnr) = F (r) . As a 



(g) result, Eq.© takes the form 



F(x)=Z-?Q j dtF 2 (t)~^- I dtF 2 (o -:-/)<•■(/) 



where 



2 K (cosht/z) 

7r 1 + e* 
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(11) 



(12) 



6 (x > 0) = 1, 6 (x < 0) = 0. The function rf> (t) decays ex- 
ponentially at t — * ±oo and has a logarithmic divergence 
at t = (see Fig. 1). For large x the last term in the 
right-hand side of Eq. lfTTj) can be neglected. After that, 
Eq. (|llj) is transformed into a differential equation which 
can be easily solved. As a result, we find the following 
solution for the screening function F: 



F(r) 



Z 
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(13) 



r>d. 

This logarithmic screening of the Coulomb potential 
results from a creation of electron-hole pairs in the vicin- 
ity of the impurity, or, in terms of quantum electrody- 
namics (QED), a "vacuum polarization" (ljj, [2oT| . This 
effect can be qualitatively described in QED by an ap- 
proach which is very similar to the Thomas- Fermi theory 
used here[20|. 

As a result, at distances much larger than the lattice 
constant charge-impurity potential in undoped graphene 
equals 



V(r) 



1 
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(14) 
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does not depend on the impurity charge Z and is much 
weaker than the bare potential Vq (r) . This follows from 
the fact that the "effective fine-structure constant" for 
graphene, a is close to 1, instead of 1/137 in QED. 

Consider now a generic case of doped graphene. In 
this case, Eqs.([T|), (O, (g]), and © result in the following 
integral equation for the total impurity potential: 



V(t)- 
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,V(v>) 
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If one neglect the nonlinear term in the right-hand side of 
Eq. (TT5]) this equation can be easily solved by the Fourier 
transform; the result for the Fourier component of the 
total potential, v (q) readsfH, 0]| 



charge impurity potential in graphene. This weakens es- 
sentially this scattering mechanism since Q In is of or- 
der of ten for typical charge carrier concentrations. Per- 
turbative estimations of the electron mobility [l2T] should 
be thus multiplied by this factor squared. As a result, 
the mobility for the same parameters turns out to be 
two order of magnitude larger. Instead of concentration- 
independent mobility, we obtain a mobility proportional 
, V 2 (r') to In 2 (kFO). This weak dependence on the charge-carrier 
concentration is probably consistent with the experimen- 
tal data[21|. More accurately, one should use an expres- 
sion for the mobility obtained by Ando[l3| (see Eq.(3.27) 
and Fig. 5 of that paper) but with the replacement of 
Z by Z* when calculating the strength of the Coulomb 
interaction. 



(15) 



where 



4e 2 M 



(17) 



is the inverse screening radius proportional to the Fermi 
wave vector kp- After inverse Fourier transformation one 
finds 



t/(r) = ~l 1 "T |HoW ^ o(Kr)1 } (18) 

with asymptotic behavior 

Ze 2 1 



V(r) 



(19) 



at /«r > 1; here Ho and Yq are Struve and Neumann 
functions. 

Estimating different terms in Eq. (|15p one can demon- 
strate that the solution ([13]) is still correct for nr < 1 and 
the solution (fl~8|) - for nr 3> 1, but with a replacement of 
Zby 



(20) 



in the latter case. Analyzing contributions to the integral 
in the right-hand side of Eq. (flU)) from these two regions 
we obtain our final result 



(21) 



1 + ZQ In 
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This is the effective charge of impurity in graphene at dis- 
tances much larger than the lattice constant. Since we 
always have kpa <C 1 this means that it is the charge that 
determines electron scattering by a long-range part of 
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